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Abstract: Super-leading logarithmic terms have previously been observed in non-global QCD 
observables. In this paper we re-calculate the first super-leading contribution to the 'gaps 
between jets' cross-section using a diagrammatic fixed order approach. We calculate for the 
first time the super-leading contribution at fifth order in a s . 



Keywords: 



^(JD,Jets 



Contents 



[l]. Introduction 



[2]. Fixed Order Calculation 



27T1 Overview and Notation 



2~2] Forming Uncut Diagrams 



2T3| Cutting and Ordering the Diagrams 



2~4\ Forming the Integrals 

2.4. 1| Super-Leading Logarithms 



Results 

Checks 



3.1 



f§ Outlook 

[S]. CH — h Implementation 

[A. 1| Forming Uncut Diagrams 

A.2| Cutting and Ordering Diagrams 



A.3| Integrals 



A. 4 Performance and Parallelisation 












U 

12 

n 
n 

H5| 



1. Introduction 

Hard QCD interactions are typically accompanied by a large number of soft and collinear gluon 
emissions. These give rise to singularities, which cancel between real and virtual emissions, 
leaving a finite remainder. However when the real emission into a certain phase space region is 
restricted, this cancellation no longer occurs, and we are left with large logarithmic corrections 
to the hard process. 

For the processes we are interested in, the leading behaviour is obtained using the limit of 
soft, wide-angle (i.e. non-collinear) gluons. We can therefore work in the Eikonal approximation. 
For virtual gluons in the Eikonal limit, we may integrate out one component of 4- momentum. 
This contour integral encloses a pole where the virtual gluon goes on-shell, and once integrated 
the remaining 3-momentum integral is of the exact form of a real gluon phase space integral, 
but with opposite sign. In an inclusive observable the cancellation between the real and virtual 



gluons in this limit is exact at the differential level. It therefore makes sense to formulate the 
calculation of this virtual gluon contribution as a phase space integral, and in particular to 
assign them rapidities corresponding to the rapidities of the real gluons they cancel against. We 
will therefore talk of virtual gluons occupying certain regions of phase space, meaning the limits 
on the remaining integrals over the virtual gluon momentum. 

An example of a restricted emission process is gaps-between-jets, where we explicitly limit 
the transverse energy emitted between hard di-jets. Traditionally [1-5] it was assumed that the 
real-virtual cancellation was perfect outside the gap so that one need only consider the virtual 
gluons lying in the gap with above the transverse energy cut-off Qq, the region in which real 
emission is vetoed and therefore cancellation is spoiled. These terms could be re-summed to all 
orders, to give a correction to the hard cross section. 

However Dasgupta and Salam [6,7] realised that one must also consider secondary emissions. 
Real gluons lying outside of the gap may emit back into the gap, and these emissions too are 
vetoed above Qq. These give rise to a new tower of logarithms, formally as important as 
the previous primary emissions, known as non-global logarithms. One must now consider any 
number of emissions, real or virtual, outside the gap regions. The colour structure for this 
quickly becomes uncalculable and hence corrections have been calculated to all orders only in 
the large N limit. 

Considering the specific case of only one gluon outside the gap, dressed with virtual gluons, 
but using the exact iV structure, Forshaw, Kyrieleis and Seymour [8] discovered a new class 
of logarithm. These logarithms appear first at order a\ relative to the hard process and are 
formally more important than the primary emission logarithm at this order, and therefore known 
as super-leading logarithms. They are however sub-leading in colour, hence were not observed 
with the original Dasgupta-Salam non-global logarithms. Their origin is traced to a breakdown 
of colour coherence in the initial state originating from the imaginary part of loop integrals, 
known as Coulomb gluons. They are therefore only present in processes with at least two 
initial-state partons and at least two final-state partons. The calculation of [8] is limited by 
the requirement to calculate increasingly large anomalous dimension matrices as the number of 
external coloured partons increases, though recent steps have been made in simplifying it [9] 
using a colour basis independent approach [10, 11]. 

In this paper we present a diagrammatic approach to calculating the gaps-between-jets 
cross-section at fixed order. This allows us to retain the exact N structure and consider any 
number of gluons outside the gap. The order that can be reached is limited by the large number 
of diagrams involved, and ultimately by the availability of computer power and time. 

In Section ^| the diagrammatic algorithm is explained. Appendix A contains some details of 
the implementation of this algorithm in C + +. In Section |^ we present results and in Section |4] 
an outlook. 



2. Fixed Order Calculation 



2.1 Overview and Notation 



The calculation uses a representation of diagrams. These diagrams may be either 'cut' or 'uncut' 
with cut diagrams representing squared matrix elements, and uncut diagrams representing a set 
of cut diagrams. Figure |I] shows how a cut diagram corresponds to the matrix element squared 
formed from a matrix element and its complex conjugate. Figure ||] shows an uncut diagram, 
representing the sum of all possible cuts. A gluon running over a cut represents the interference 
of two real gluon emission diagrams. A gluon to one side of a cut represents a virtual gluon. 
All cuts through a diagram share the same colour factor (up to a possible overall sign, which 
we discuss in detail later). 

We work in the Eikonal (soft) gluon approximation, in which all components of the gluon 
momenta fcj are much less than any components of the momenta Pi of the original hard partons. 
In this limit, one need only consider Feynman diagrams in which the soft gluons are connected 
to the external hard partons, and not to internal lines within the hard scattering. Replacements 
of the form 

yy+jO ^ p m (2 1} 

(p + k) 2 + ie p.k + ie 
are made in the expression for the matrix element. Equation Q2.1| ) represents the emission of 
a gluon of momentum k from a quark of momentum p, however the form of the right-hand 
side is universal whether the parton is a quark or gluon. In the squared matrix element for 
gluon emission, where there are Eikonal attachments at both ends of the gluon (or in the matrix 
element for virtual gluons where this is always true), all gluons give terms of the form (in 
Feynman gauge) 

(2.2) 



Pi-Pi 



(pi.k ± ie)(pj.k ± ie) 

where pi and pj are the momenta of two partons joined by a gluon of momentum k. The ± 
in the denominator depends on whether the gluon is emitted or absorbed, and whether the 
emitters sit in the initial or final state. For triple gluon vertices a further complication exists 
due to the anti-symmetric nature of the colour- and momentum-dependent parts of the triple 
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Figure 1: A cut diagram represents a squared matrix element. 



gluon coupling. In order to correspond to our convention (the right hand side of eq. ( |2.1| )) the 
soft gluon must be represented as being emitted to the right as one looks in the direction of 
momentum of a harder emitting gluon. 

We retain terms leading in the logarithm of Q, the energy exchanged in the hard process, 
over the cut-off Qo] the terms are also leading in Y, the rapidity separation, or 7r 2 , which arises 
from Coulomb gluon terms. This is referred to in [12] as the tc 2 Double Leading Logarithm 
Approximation (tt 2 DLLA). The colour basis independent approach of [9] shows that the real 
part of sub-leading terms in Y in the anomalous dimension matrix are proportional to the 
identity matrix, and that this approximation contains all the interesting physics, in particular 
the super-leading logarithmic terms that interest us. 

The leading logarithms are obtained using the strong ordering condition. This states that 
each emission is much softer, in some chosen kinematic variable, e.g. kx, than the last and 
can therefore be treated as independent. Each emission sees the result of previous emissions 
as the hard process and can therefore be calculated knowing only the momenta and colour 
configuration of the external partons. Previously added virtual gluons simply alter the colour 
structure of the hard process, real gluons mean the new emission is added, e.g. transforming a 
hard 2 —>■ 2 process to a 2 —>■ 3 process. The addition of each soft gluon to an external hard 
parton leads to a logarithmic integral over the internal hard parton propagator. The nesting of 
multiple emission results in a leading contribution with one logarithm for each order of a s . 

Colour coherence would mean that the choice of ordering variable is unimportant: strong 
ordering in any kinematic variable of the soft gluons would give the same result. However as 
the appearance of super-leading logarithms is interpreted as a breakdown of colour coherence, 
the variable used may be crucial. In common with previous calculations, we use the transverse 
momentum as our ordering variable (Q 3> &t,i 3> ^t,2 3> ■■• 3> &t> ^qcd)- This is, however, 
an area of on-going investigation. 

It should be noted that, in the various approximations we use, all cuts through the softest 
gluon, as in figure |2|, cancel. When a phase space restriction is introduced this cancellation fails. 
This only applies for the softest gluon - when there is more than one gluon, complicated colour 
structures can ruin this cancellation, and one must consider all cuts through all diagrams. 

By including all topologically distinct uncut diagrams, which include all allowed cuts, we 
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Figure 2: An uncut diagram represents a sum of cut diagrams. 



Figure 3: The set of all uncut first order diagrams. 



can build up the basis of the entire order as correction to the hard process. 



2.2 Forming Uncut Diagrams 

The calculation starts with an original hard process (e.g. qq — > qq exchanging a gluon). The 
hard process is formed into an uncut diagram, that is, combined with its complex conjugate to 
form a closed loop. Only when this diagram is 'cut', i.e. separated into a matrix element and 
complex conjugate pair can it be calculated physically. 

From these uncut diagrams higher orders are built by iteratively adding gluons attaching 
to all possible external legs and previously added gluons. Real gluons may re-emit a gluon, 
although coupling to virtual gluons (including any in the initial hard process) is suppressed. In 
the Eikonal approximation using Feynman gauge, gluons attached at both ends to the same leg 
make no contribution, so these are not considered. 

The colour factors corresponding to these diagrams are calculated at this stage. Here we 
must choose a representation of the triple gluon vertex, which may introduce a factor of —1 once 
the momentum direction of the partons is set later on. This process of adding gluons leaves us 
with a set of uncut diagrams each with a colour factor, a polynomial in the number of colours 
N. Figure [5] shows all the possible first order uncut diagrams attaching to a qq — > qq process. 



2.3 Cutting and Ordering the Diagrams 

The next stage is to cut these diagrams. All possible cuts are considered. All original hard lines 
must be cut once and only once, leaving a hard scatter either side of the cut. Vertices including 
added gluons, both those attached to hard partons and triple gluon vertices may sit either side 
of the cut. 

The cut diagrams are subsequently ordered. We assume strong ordering, in order to capture 



the leading logarithms, as explained in section [2.1| . All restrictions are obtained by moving 
outwards along the hard parton lines from the hard process to the cut, and similarly along 
added gluon lines. Diagrams with more than one possible ordering have all orderings retained 
separately, those which are inconsistent with strong ordering are rejected. Also rejected are cuts 
leading to disconnected gluons and those where soft gluons attach only to one harder parton (it 
can be seen from eq. (|2.2|) that these terms will be for massless partons). 



Figure 4: All possible cuts through a second order diagram. 



Figure | shows all possible cuts through a second order diagram. Note that diagrams d.) 
and f.) are the same, but two possible orderings exist so we keep a diagram for each. We can see 
in diagrams a.) to c.) the leftmost gluon must be softest, as it is furthest from the hard process 
to the right of the cut. We might also notice that diagrams a.) to d.) resemble figure 0, and 
we are free to assume the leftmost gluon can be softest in diagram d.) as well. This set forms 
a cancelling set of all cuts through the softest gluon. In order to include all possible orderings 
diagram f.) must be the alternative ordering of d.) with the rightmost gluon being softest, and 
forming a canceling set with g.) to L). To the left of the cut in e.) the leftmost gluon is closest 
to the hard process and therefore hardest. But to the right, the rightmost gluon is hardest. 
This arrangement is incompatible with the strong ordering we assume. Diagrams e.) and j.) 
must therefore be rejected. 

At this stage triple gluon vertices are also considered - the emitting gluon must be real, and 
harder than the emitted gluon and we impose these conditions during the ordering phase. Three 
gluons meeting at a vertex, treated equally when obtaining unique uncut diagrams now become 
one hard gluon emitting one softer gluon. The ordering stage determines the momenta of the 
gluons, and we may now resolve the factor of —1 ambiguity resulting from making a specific 
choice of order in the triple gluon coupling, ifijk when we calculated the colour factors of the 
uncut diagrams. 

For example in figure [|, the colour factor is calculated at the uncut diagram stage, shown on 




Figure 5: A triple gluon diagram and two possible orderings arising from one cut. 



the left, representing the triple gluon vertex as 1/123. The first cut diagram shows one possible 
ordering with gluon 1-2 emitting gluon 3. Momentum flows from the emission of gluon 1, and 
gluon 3 is emitted to the right as one looks along the momentum direction, the convention we 
have adopted for the Eikonal approximation. The second cut diagram shows an alternative 
ordering of the same cut, with gluon 3-2 being hard and gluon 1 soft. Momentum flows from 
the emission of gluon 3 upward, and gluon 2 is emitted to the left of the momentum direction. 
As our convention requires the colour factor to be written as if the gluon is emitted to the right, 
in this case, if 213, we acquire an overall minus sign for this second cut diagram, which is taken 
into the colour factor. 



2.4 Forming the Integrals 

Now we have all the information needed to perform the calculation, a set of integrals over the 
added gluon four momenta is formed from these ordered diagrams. Each integral is defined only 
by the gluon whose momentum we integrate over, the two partons the gluon stretches between 
(of the form in eq. (2T2")) and whether the gluon is real or virtual. For example, if we take the 



cut diagram b.) from figure ^ we have a harder virtual gluon, k\, and a softer real gluon, k 2 , 
and from this we form: 

°"o ^, An 4 f V3-Va i f p 3 .p 4 d 3 k 2 



° = TT7TrMN)gl 



C {N) '* s J virt (ps.h + ze)(p 4 .fci - ie) k\ + le (2vr) 4 J real (p 3 .k 2 )(pA-k 2 ) 2fc°(27r) 3 

(2.3) 

where a and C are the cross section and colour factor of the original hard process, C the colour 
factor of the new process and p 3 and p^ the momenta of the two outgoing quarks. 
Writing the momenta in the centre of mass frame of the incoming partons, we get 



Pi 


= 0,0,1) 




(2.4) 


P'2 


= ^(l; 0,0,-1) 




(2.5) 


P3 


= Q (cosh(cfy/2); 0, 1, sinh(<5r//2)) 


(2.6) 


Pi 


= Q(cosh(<fy/2);0,-1, 


-sinh(&7/2)) 


(2.7) 


h 


= k T ,i (cosh(yi); sin(0i) 


,cos(0i),sinh(yi)) 


(2.8) 


k 2 


= k Tt2 (cosh(2/ 2 );sin(0 2 ) 


,cos(0 2 ),sinh(2/ 2 )) 


(2.9) 



The dk\ integral in the virtual gluon can be evaluated as a complex integral and gives two poles, 
one corresponding to gluon 1 being on-shell (known as the eikonal gluon), and one corresponding 
to an emitting parton being on-shell (known as the Coulomb gluon) giving an imaginary term in 
if the emitters are both in the initial or both in the final state (as they are in this example). After 
performing the azimuthal (0) integrals it can be shown that the rapidity dependent integrands 




Figure 6: An rj-cj) plot showing the jets (red circles) and gap (grey band). 



are of the form 



2sinh(|j/|) 



and in the large \y\ limit reduce* to a step function 9(y). This leaves 



eq. 



as 
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dy 2 



(2.10) 



where Qo is the energy cut-off into the gap, kr,2 is cut off by as we have assumed strong kx 
ordering and the real gluon rapidity is restricted to the region outside the gap (see figure H). In 
general these contributions come with an overall sign that must be determined from the colour 
factor, the number of cut gluons, the orientation of triple-gluon vertices, etc. 

As we can see, the rapidity integrals only depend on the limits. We also keep Coulomb 
gluon terms, which are the imaginary part of the virtual gluon integrals of the form Z7T. The 
physical observable is real, but may contain terms proportional to n n arising from these terms, 
where n is even. 

We notice only the colour factor and rapidity integrals change between diagrams, the trans- 
verse momentum and numerical factors change only between orders in a s . After azimuthal 
averaging in the high energy (large Y) limit the rapidity of the gluon is restricted to lie between 
that of the two partons it connects. Therefore, for each gluon, we get a rapidity dependent 
contribution dependent only on which region on the plot ^ it occupies. Gluons coupled to other 
gluons form a set of nested integrals. Answers are given in terms of these rapidity zones, the 
width of the gap region, Y, the out-of-gap region, Y, and the jet radius, R = \{5rj — Y). 

Out of gap contributions (Y terms) cancel over the set of all diagrams up to order cuf. When 
this cancellation fails, the out of gap region (with the rapidity of the incoming partons infinite) 
must be dealt with specially, as described in the next section. 

*Thcre is a subtlety here: the original expression is divergent at the point y = 0, corresponding to emission 
collinear with one of the hard partons, while the approximated expression is finite (although discontinuous) 
there. The approximation is justified by the fact, proved using the colour basis independent notation in [9], that 
final-state collinear divergences cancel to all orders and so do not have an effect on the final result. The large Y 
limit is therefore captured by this step function approximation. 



2.4.1 Super-Leading Logarithms 

The out-of-gap rapidity region appears to be infinite, but at large rapidities a collinear rather 
than soft approximation is more appropriate. Matching the soft and collinear regions in [8] it 
was argued that the out-of-gap rapidity interval should be cut off at Y — 2 In ■ These are 
the super- leading logarithms, which must slot into the nested kx integrals (as seen in eq. ( 2.1(J| )). 

We require an extra factor, depending on where in the k? ordered chain the gluon lies. 
With no out-of-gap contribution we get something of the form ln n (j^j /n\. For each additional 
logarithm from an out-of-gap gluon, n increases by one in both the power of the logarithm and 
the factorial, and there is an extra numerical factor of the number of harder logarithms plus 
one. For one gluon outside the gap at fourth order, for example, we have 



21n 5 f# 

\Qo 



where m is the number of gluons harder than the out-of-gap ( "super- leading" ) gluon. Similar 
equations hold for larger numbers of out-of-gap gluons. 

We can now re-create the full answer by taking the summed colour and rapidity dependent 
parts calculated on a diagram-by-diagram basis, and combining them with the kx logarithms 
(accounting for superleading logarithms), as and hard cross-section pieces calculated on an 
order-by-order basis. 

We have written a C + + program to encode the calculation, in order to systematically 
ensure all diagrams and cuts were accounted for. Further details are contained in the appendix. 



3. Results 

Here we present the results as corrections to <jq, the hard scattering process with no additional 
soft gluons. The main original results are the two gluon outside the gap coefficients (eqs. ( |3.5| ) . 
( p,10| ), ( |3.15| ), ( |3.20| )), but we show the full it 2 DLL A results at fourth and fifth order in a s , 



calculated at the same time for comparison. The one gluon outside the gap coefficients at fifth 
order, and for colour- singlet hard processes, are also new. 

We present here results for 0,1 and 2 gluons outside the gap. The calculation includes any 
number of gluons outside the gap, contributions from higher numbers of gluons outside cancel 
in the sum over diagrams. This explicitly shows, for example, that there is no two gluon outside 
contribution at O(a^), cis £ir erued in [81. 

The results given are for the ^-channel scattering of various particles exchanging a gluon, 
and for two quarks exchanging a colour singlet. Although the colour singlet case is sub leading 
compared to the gluon exchange for simple scattering, it shares a colour structure with Higgs 
production through vector boson fusion, for which a rapidity gap is an important experimental 



signal, and which has been the subject of previous studies into soft gluon resummation [13]. 
Exchanging an anti-quark for a quark leaves the results unchanged. 
For quark-quark scattering, exchanging a gluon at order a 4 : 
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For quark-gluon scattering we find: 
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For gluon-gluon scattering we obtain: 
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For quark-quark scattering exchanging a singlet: 



,"2aV 4 /Q\ 2 ^ 2 iV 2 -l 
<Wj = a [—J ]n \^-JnY -32-» ( 3 - 16 ) 

<W< = ( — Y'l'."' (-^| vr 2 y^-l, (3.17) 



/2a s \ 5 5 / Q 
<ro, qsq , a s = -ao ^— j In 



7r y 7r y 

960 240iV 



(3.18) 



It is interesting to note that in each case the fifth-order one-outside and two-outside cross 
sections have the same sign, opposite to that of the fourth-order one-outside cross section. 

3.1 Checks 

A series of checks against pre-existing results were undertaken in order to validate the perfor- 
mance of the algorithm, and gain confidence in the original results. 

The number of unique uncut diagrams produced is in theory calculable, however it is very 
complex for the large numbers of diagrams we dealt with. However by isolating subsets of 
diagrams (e.g. for two or three hard parton states, for five gluon configurations containing 
seven triple gluon vertices (and therefore only three hard attachments) over four hard partons) 
and comparing with the qgraf program [14] we could be confident we had the full set of unique 
diagrams at each order. 

A major complication encountered in the completion of this work was the correct assignment 
of the ±1 factors arising in triple gluon vertices. The failure of this aspect of the calculation 
gave rise to imaginary (in) terms remaining in the final answers and we consider the eventual 
elimination of these terms a good test of our sign convention. 

Terms with two or more gluons outside the gap, although not explicitly calculated previously, 
have been argued to cancel below 0(a 5 s ) [8] and our results show this explicitly. We also find 
those with three or more gluons cancel at 0(«f) as expected. 



The leading Y and tt terms reproduce the result of a ladder of virtual gluons lying in the 
gap and have been checked against other calculations at all orders. The one gluon outside the 
gap terms reproduce the results presented at order aj in [9], and have been compared for the 
quark-only processes at order a b s with a further expansion of the exponential approach of [8]. 

4. Outlook 

In this paper we have presented a diagrammatic method for calculating large logarithmic cor- 
rections arising from the addition of soft gluons to a fixed order in a s . We have presented the 
first two gluon out of the gap corrections to the gaps-between-jets cross-section. The results 
support the speculation of [8], that at each higher order one obtains an additional logarithm 
from contributions with an additional gluon outside the gap, and hence that the true leading 
behaviour is not a™ log™ Q/Qo, as once thought, but a"log 2n ~ 3 Q/Qo- Clearly to calculate such 
a contribution to arbitrary orders, one needs control over the colour structure of arbitrarily 
complicated final states. 

As previously mentioned, one of the on-going areas of investigation in order to understand 
the origin of the super-leading logarithm terms is a study of the choice of ordering variable. 
Although we present here results obtained using transverse momentum, our approach only 
requires this to be specified at a late stage of the calculation. The full set of ordered diagrams 
could be used with a series of integrals calculated in a different ordering scheme in order to test 
the assumptions made. 

The presence of super-leading terms in other non-global observables is expected and, as for 
the ordering variable, the ordered diagram sets remain valid. To use these one must relate the 
structure of the diagrams to the integrals for the observable required. 
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A. C-\ — |- Implementation 

A.l Forming Uncut Diagrams 

The unique set of uncut diagrams was achieved by iteratively adding a gluon in all possible ways 
to the previous order set of diagrams. To ensure only topologically distinct sets were retained 
we used the C+ + standard template library set template. For this we needed to define a unique 
numbering system for each diagram, and a less-than operator, which allowed the set functions 
to identify unique diagrams. 



The gluons were numbered as follows. Gluons attached to one of the hard partons were 
numbered first, starting at one, increasing from left to right along the hard lines in turn. Then 
internal gluons are dealt with, the lowest numbers remaining are assigned to the gluons joining a 
vertex containing the lowest already numbered gluon, e.g. if n gluons attached to hard partons, 
two gluons meeting gluon 1 at a triple gluon vertex would be assigned the numbers n + 1 and 
n + 2. To decide which of the two gluons receives which number, the lowest number is assigned 
to the gluon closest to gluon 1 (i.e. with fewest vertices in between). If the two gluons are 
equidistant from gluon 1 (as in our previous example) their distances to gluon 2 are compared 
and so on until they can be separated. This continues until all gluons are numbered. 



Figure 7: A fifth order diagram numbered according to our l° west numbers attached, 4 on the 



line. After that the vertex joining 2 is done, with 6 numbering the line closest to 1. Next we 
treat the vertex containing 3, gluon 7 is already assigned so the remaining gluon is labelled 8. 
We then apply the same algorithm to with the vertices containing 4, 5, 6 and all gluons are 
numbered. The diagram in figure would be stored in the program as: 



Line 0: 1 




Line 1: 2 


3 


Vertex 0: 


14 5 


Vertex 1: 


2 6 7 


Vertex 2: 


3 7 8 


Vertex 3: 


4 9 10 


Vertex 4: 


5 10 11 


Vertex 5: 


6 9 12 


Vertex 6: 


8 11 12 



The number of diagrams we would expect, corresponding to our constraints, could not be 
trivially calculated. However we compared numbers for subsets of diagrams (e.g. two and three 
parton states) using the qgraf program [14] and found agreement in all the areas we checked. 
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We now have a unique descriptor 
of each diagram consisting of a list 
of lines with numbered gluon connec- 
tions, and a list of triple gluon ver- 
tices, each with three attached glu- 
ons. An example is shown in figure |7]. 
First gluons 1, 2 and 3 are numbered, 
as they join hard lines. Then the ver- 
tex containing gluon 1 has the next 



scheme. 



line closest to 2 (three vertices sep- 
arate the two), 5 on the remaining 



At this stage in the code, the colour factor of each uncut diagram was calculated via an 
interface to the COLOUR program [15]. For the triple gluon vertices we simply entered the 
connections in increasing numerical order, as stored by the code's internal numbering system. 
The factors of —1 originating in the anti-symmetric nature of the triple gluon coupling were 
resolved later, as described in section |273|. The colour factors are stored as polynomials in N 



A. 2 Cutting and Ordering Diagrams 

The full set of cut and ordered diagrams was obtained by systematically producing all cuts and 
orders and discounting those that failed our requirements. Connections on hard lines and triple 
gluon vertices were simply assigned in turn to either the left or right of the cut. Diagrams with 
disconnected gluons were thrown away at this stage. 

The ordering constraints were kept track of through a series of theta functions, each describ- 
ing a less than-greater than relation. Moving outwards from the hard process along the hard 
parton lines strong ordering was recorded, when diagrams produced contradictions they were 
thrown away. Then for remaining diagrams the triple gluon vertices were treated in turn, with 
one gluon chosen as the softer gluon, and the remaining two joining to form the hard emitted 
gluon. Resulting constraints were added at each stage and diagrams discarded as soon as any 
inconsistencies were found. 

When all triple gluon vertices had been dealt with, the result of the fully ordered diagram 
was stored as a list of gluons with information on which hard partons or other gluons they 
connected, and whether they were virtual or real. 

A. 3 Integrals 

In the high energy limit, there are very few required integrals, and they were calculated by hand 
(see section p.4j) and only the colour factors and rapidity integrals change between diagrams. 
The definition of the observable we calculate is most easily described in terms of a cone algorithm 
with radius R, but in fact since we use an approximation in which all emissions are assumed to 
be well separated in rapidity our results are general for any jet algorithm. 

The main complication that exists is that the rapidity integral limits when there are many 
triple gluon connections lead to a series of complex nested integrals, which are dealt with as 
follows. 

A series of theta (step) functions describing the rapidity limits are taken, describing which 
of the rapidity regions illustrated in figure | the gluons fall into. 
So, for example, from equation ( j2.10| ) we would retain 



C(N) 



o[ S 4-m)o( yi + S 4 



2 a V 2 



ei S 4-»*)oL + s 4 



2 a V 2 



(A.1) 

y2=real 



yi=virt 

plus an imaginary part. Terms not noted here are common to all diagrams at this order. 



The virtual gluon produces a term Y + 2R (where R = Sr) 2 Y is an effective jet radius) and 
the real gluon a term 2R (being vetoed in region Y). So the total rapidity dependent part is 
2RY + AR 2 . 

If a gluon lies between, for example, —Sr]/2 and Y/2, it is split into two terms, one with an 
in-gap gluon and one with an out-of-gap gluon 

and the resulting terms are Y + 2R and Y . 

Gluons may also be re-emitted, and therefore limited by the rapidity of the emitting gluon, 
which produces nested integrals. For example 

6 - vi) 6 + vi) 9 - y*) 6 d + v*) %i - yi) ( A - 3 ) 



s 2 J \2 J \2 J \2 
gives Y 2 /2\. To deal with all cases of this kind, we expand the theta functions, so for example 

- j/i Wtt + i/i - y*) d (^ + y2 



k 2 J \2 J \2 °'J \2 
^e(j-y i y(^+y i y(^-y 2 y(^+y 2 S \ 9{ Vl - y 2 ) 

+6 (| - yi) (| + y?j 6 (| - y^j 6 (| + y^j 6(y 2 - Vl ) (A.4) 

giving a result of Y 2 /2 + Y 2 /2 = Y 2 as we would expect for two gluons spanning the gap region. 
However this expansion has the advantage of automatically including all term of the form of 
eq. (TAJ). 

With these fully expanded set of theta functions the rapidity dependent integral is the width 
of a rapidity region raised to the power of the number of gluons in that region, over the factorial 
of that number of gluons, summed over all rapidity regions: 

4^ rn\ 

i 

where (Ay)j are the width of the regions, i, (e.g. Y, R, Y) and rij is the number of gluons lying 
within each region. 

This simple formula holds for all diagrams after the theta function expansion has taken 
place. The performance suffers due to the large increase in the number of terms to be dealt 
with, but the benefits are the simplicity of the integration, and being sure we have systematically 
included all the required terms. 

A.4 Performance and Parallelisation 

The performance of the code decreases rapidly as the 0(a s ), and therefore the number of 
diagrams, increases. The table below shows the number of diagrams at three stages of the 
program, and the approximate running time (including compilation time). 



0(a s ) 


Uncut diagrams 


Cut diagrams 


Ordered diagrams 


Run time(s) 


1 


6 


24 


24 


21 


2 


55 


564 


672 


22 


3 


681 


17,194 


30,048 


44 


4 


10,529 


529,910 


1,746,272 


3700 


5 


193,741 


22,521,544 


127,424,961 


1,600,000 



It can be seen that the processing time increases at a higher rate than the number of uncut 
diagrams, as the number of cut diagrams, possible orderings of those diagrams and (not shown) 
expansions to account for nested integrals (section |A.3|) all increase per-diagram with order. 

Although times were manageable up to fourth order, to reach fifth order it was necessary to 
run the program in parallel sessions. As the uncut diagrams need to be checked to ensure they 
form a unique set, they must be created for each order as one stage. The uncut connections 
were then written to a series of text files. The connections were constant for any 2^2 process, 
so this stage only need be done once. 

Each file was read in and the pre-prepared connections combined with a 0th order diagram 
(qq — > qq, gg — > gg etc) and colour factors calculated. From here we have the normal set of 
uncut diagrams and the code can proceed as usual with the following steps. Results in the form 
of a colour factor (polynomial in N) and term in Y, Y, m and R are then written to a text file. 
The program can then read in these results files, add them and output the result. 

A similar process can be used within each run of the program. The uncut diagram set is 
split into smaller sets, and these sets run through the intermediate stages to a result. A small 
number of uncut diagrams can quickly give rise to a large number of cut diagrams, ordered 
diagrams and integral terms, which can cause problems in temporary memory, before collapsing 
down to a small number of result terms. Subdividing allows the full result for the larger set of 
uncut diagrams to be calculated without at any one time storing a large number of the diagrams 
necessary for the intermediate stages. 
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